Formation energy in a-phase Fe-V alloys 
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Formation energy of the cr-phase in the Fe-V alloy system, AE, was computed in the full compo- 
sitional range of its occurrence (~ 34 < x <~ 60) using the electronic band structure calculations 
by means of the KKR method. A_E-values were found to strongly depend on the Fe concentration, 
also its variation with different site occupancies was characteristic of a given lattice site. Calculated 
magnetic, Smagn, and configuration, S CO nf, entropy contributions were used to determine sublattice 
occupancies for various compositions and temperatures. The results agree well with those obtained 
from neutron diffraction measurements. 

PACS numbers: 71.15.Nc 61.43.-j, 71. 20. Be, 71.23.-k, 74.20.Pq, 75.50.Bb, 
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I. INTRODUCTION 

Sigma (cr) phase was found in over 40 binary alloys and 
many other three- or multi-component alloy systems^. 
The phase cannot be formed at the stage of solidifica- 
tion of the solution of alloying elements. Instead, it can 
only be obtained by a high temperature annealing pro- 
cess (solid state reaction). From the viewpoint of tech- 
nological applications of alloys, the cr-phase should be 
avoided as its presence drastically deteriorates different 
mechanical properties of materials in which it precipi- 
tates. A better knowledge of conditions of its formation 
may be helpful for fabrications of materials in which the 
unwanted phase does not precipitate or its precipitatioin 
is significantly retarded. On the other hand, due to its 
complexity (30 atoms per unit cell, five different sublat- 
tices, high (12-15) coordination numbers, and lack of a 
stoichiometry) make the cr-phase a very challending ob- 
ject for theoretical calculations. They are very helpful for 
a proper interpretation of measurements, especially those 
done by means of microscopic methods like Mssbauer 
spectroscopy or nuclear magnetic resonance, as well as 
for a better understanding of the results obtained^— 

In the Fe-V system the cr-phase can be formed by an 
isothermal annealing of the free-master alloys in a wide 
range of concentration (~ 34 — 60 at.% V) and tempera- 
ture (~ 400 — 1230°C)i. Once formed, it 'remains stable 
at lower temperatures, but it can be dissolved into the a- 
phase at temperatures above ~ 1230°C. Mechanism of its 
formation as well as that of its dissolution is not known 
with sufficient precision yet. Knowledge of the energy of 
its formation compared to the energy of formation of the 
ct-phase (from which it precipitates) may help to better 
understand processes responsible for its formation and 
dissolution. Thus, it may be useful in a fabrication of a 
new generation of materials, like stainless steels, having 
suitable properties for a construction of new generations 
of industrially important facilities. 

Having a tetragonal unit cell and hosting 30 atoms 
distributed over five crystallographically non-equivalent 
sites (usually termed A through E) this phase remains 
still challenging for theoretical calculations. Luckily, in- 



crease of computing capabilities of modern computers 
makes it possible that more sophisticated models can be 
applied to more adequately and precisely describe the 
complex structure and properties of the cr-phase^—. 



II. COMPUTATIONAL DETAILS 

The present analysis of the structural stability of 
the cr-phase was based on electronic structure cal- 
culations, similar to the methodology applied for cr- 
FeV and cr-FeCr alloy systems to study hyperfine 
interactions^&ii as well as the formation energy for the 
latter—. We have employed the Korringa-Kohn-Rostoker 
(KKR) technique^— to calculate the electronic struc- 
ture and the total energy in the framework of the non- 
relativistic LDA theory and using the muffin-tin approx- 
imation for a crystal potential. More computational de- 
tails on cr— phase can be found in Refs. d and|j. 

Since in the real structure of cr-FeV all five crystal- 
lographic sites are occupied by both types of atoms, a 
chemical disorder on these sublattices must be taken into 
account. A more sophisticated approach based on the 
KKR method with the coherent potential approximation 
(CPA), allowing to treat the disorder as random, should 
be capable of accounting for 2 30 atomic configurations, 
a target that still remains unreachable. In these circum- 
stances, an alternative, but a feasible approach consists 
in selection of some ordered elementary unit cells with 
different configurations of atoms on the five sublattices, 
which permits to use the KKR method (instead of KKR- 
CPA one). In practice, the electronic structure calcula- 
tions were performed for a finite number of the unit cells 
with a symmetry reduced to a simple tetragonal one (PI 
instead of P42/mnm), in which each position (but not the 
sublattice) was occupied by one type of atoms (Fe or V) . 
The atoms were distributed over the five different sites 
with the probabilities determined from neutron diffrac- 
tion measurements^ which were next altered in the way 
described below. Differences between total energies (per 
atom), AE (usually called the formation energy), deter- 
mined for each of these atomic configurations, E a , and 
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corresponding values obtained for pure constituents, i.e. 
bcc-Fe and bcc-V, E a , could be calculated^ according to 
Eq. ©: 

AE = E a -{\- x)E^ - xE v a (1) 

where x stands for the concentration of V. They can be 
further presented in form of a diagram as a function of 
concentration of the alloying elements. This approach 
was already successfully applied to the analysis of the for- 
mation energy of the cr-phase in the Fe-Cr alloy system 5 . 

As mentioned above, the cr-phase is characterized by a 
chemical disorder on all five sublattices, and the concen- 
trations of Fe/V atoms on each of them were determined 
experimentally- 5 ". Since both Fe and V atoms are present 
on all five sublattices, it is useful for further considera- 
tions to introduce shortened notations for samples with 
different sublattice occupancies. Multiplicity (denoted 
by Ni) of the five Wyckoff positions, namely 2a, 4f, 8i, 
8i' and 8j correspond to A, B, C, D and E site, respec- 
tively. Since the total number of atoms on each of the 
five sublattices was already known, thus it was enough 
to indicate the number of the Fe atoms on each sublat- 
tice explicitly, while the V atoms constitute the balance. 
According to this convention, an illustrative example of 
cr-20172 describes a system comprising 2, 1, 7, and 2 
atoms of Fe on the sites A, C, D and E, respectively, and 
also 4, 7, 1 and 6 atoms of V atoms on the sites B, C, D 
and E. 

In the present work we analyzed more carefully ordered 
unit cells corresponding to the following atomic configu- 
rations: cr-20172, cr-20383 and cr-21485. The total num- 
ber of Fe atoms in them is 12, 16 and 20, respectively. 
The three configurations have been chosen intentionally 
to best approximate experimentally investigated samples 
viz. Fe4iV59, Fe52V4s and FS66V34, respectively^. These 
three cases correspond to the border and middle com- 
positions of the cr-phase existence on the phase diagram 
of the Fe-V alloy system—. Actually, we treated these 
configurations as the reference for a further analysis of 
AE. 

In the next step we have calculated the AE- values, for 
the unit cells changed (with respect to the three reference 
ones) by a replacement of one or two atoms exclusively 
on the same site. Ten ordered approximants of each con- 
figuration of sublattice occupancies with different atomic 
arrangements were analyzed, and average AE- values, ob- 
tained in this approach for neighborhoods of the three 
reference configurations, are presented as a function of 
Fe atoms number per unit cell, riFe (Fig. [JJ). Simulta- 
neously, the Ai?-changes caused by the replacements of 
atoms on each sublattice are indicated with various col- 
ors. 



III. RESULTS AND DISCUSSION 

The changes of AE in the vicinity of the three refer- 
ence configurations are fairly linear functions of riFe fo r 
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FIG. 1: (Online color) Calculated formation energy differ- 
ences for the Fe-V cr-phase, AE, for the chosen three basic 
atomic configurations and its neighborhood, versus the num- 
ber of Fe atoms per unit cell, riFe (V concentration is also 
indicated). Solid lines stand for the best fit to the data - 
see Appendix. A_B-values determined for occupation changes 
(with respect to the reference configurations) on the five lat- 
tice sites are indicated by different colors. 



the sites A, B, C, and D, and rather non-linear for the 
site E. The slopes of the lines for the sites A, B and C are 
characteristic of the site and equal -0.021(1), 0.025(1) and 
0.019(1) - in eV/atom - respectively. Interestingly, these 
slopes are more or less the same for all reference config- 
urations (see, Fig. [TJ, while for the site D the value of 
the slope depends on the alloy basic configuration around 
which it is calculated (changes from -0.022 for cr-20172 
to -0.009 eV/atom for cr-21485). AE corresponding to 
the changes on the site E departs from linearity and can 
be approximated using a second-order polynomial of the 
Fe-concentration on this site. 

The model used for the formation energy calculation 
allowed to determine the values of AE only for such cases, 
in which individual sites (but not sublattices) are occu- 
pied by just one type of atoms. This really corresponds 
to a one selected unit cell, whereas in the neighboring 
cells atoms may be arranged in a somewhat different way 
on the sublattices. For that reason, for a more realistic 
description of crystalline samples it is more convenient 
to operate with the average formation energy that corre- 
sponds to the average occupancy of the sublattices by the 
alloying elements. As shown in Fig. [TJ AE have been de- 
termined for almost 50 different configurations of Fe and 
V atoms, including - to some extend - the chemical dis- 
order on the sublattices. The relatively simple character 
of the Ai?-changes described above (mostly linear behav- 
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ior), provides a possibility of determining a dependence 
of AE on the Fe-concentrations, x il i =A through E, 
on each sublattice. It can be achieved by approximation 
of AE using a function depending on the variables Xi, 
on five individual sites, AE = AE(xa,xb,xc,xd,xe)- 
Taking into account the established character of the AE- 
changes, there were assumed to be linear (A trough D 
sites), and quadratic (E site) dependences on Xi (see Ap- 
pendix). Fitting such-defined functions to the calculated 
AF-values (using the least-square method) yielded nu- 
merical values of polynomial coefficients. Consequently, 
it allowed to determine AE for any atomic configuration 
for all five sites in the a-FeV structure. 



A. Entropy contributions 

Bearing in mind that the aforementioned total energy 
E and the electronic structure analysis correspond, in 
principle, to ground state properties, one can attempt to 
discuss crystal stability at finite temperature T, applying 
commonly used scheme based on a computation of the 
free energy, F = E — TS, where S stays for an entropy. 
Considering the process of the er-phase formation, one 
can discuss AF = AE — TS, being the difference between 
the free energy of cr-FeV alloy and free energies of its 
constituents (namely, weighted contributions of a pure 
bcc Fe and V which are already known). 

The total entropy that actually should be taken into 
account in the present case has several contributions, 
namely: magnetic entropy, S mag , configuration entropy, 
Sconf, electronic entropy, S e i, and phonon entropy, S p h- 
The first two terms could be computed, and the results of 
these computations are shown in Fig. 2, for the magnetic 
entropy, and in Fig. 3, for the configuration one. 

Using the magnetic moments maximum saturation en- 
tropy for each atom on the sublattice, S™%g = kg hi(2 J+ 
1), can be determined 18 . In this work SVnag-values were 
calculated according to the following Eq. ©: 
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k B ]T(P, ln(/xf e + !) + (!- Pi) ln( M r + 1)) (2) 



where k B is the Boltzmann constant, Pi stands for the 
occupancy of Fe on each sublattice, \ii is the average mag- 
netic moment of Fe/V atoms belonging to this sublattice 
and summation (i-index) runs over all five sublattices 
i.e. from A to E. /^-values for each sublattice and differ- 
ent Xi configurations were obtained from the electronic 
structure KKR calculations. Fig. 2 gives a clear evi- 
dence that for each sample S mag increases with riFe for 
each lattice site. The increase of the magnetic entropy 
remains in line with our previous calculations reporting 
that the magnetic moments (also the hyperfine fields) 
at particular sites increase with npj±. Similarly to the 
AE expression, S mag can be represented as a function of 
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FIG. 2: (Online color) The magnetic entropy, Smag, calcu- 
lated versus the number of Fe atoms per unit cell, np e , for 
different atomic configurations (samples) of the cr-FeV. The 
concentration of vanadium is indicated, too. Solid lines stand 
for the best fit to the data - see Appendix. S mag -values de- 
termined for occupation changes on the individual lattice site 
are indicated by different colors. 



linearity. To approximate the S mag dependence on Xi, 
we assumed a mathematical form of the function analo- 
gous to that used for AE (see Appendix) . The numerical 
values of the polynomial coefficients were determined by 
fitting the function to the calculated S^^-values. 

The configuration entropy, shown in Fig. 3, was calcu- 
lated according to Eq. ([3]): 

Sconf = -kg ^(P, HP) + (1 - Pi) ln(l - Pi)) (3) 

It shows a strong dependence on the site for each of 
three basic configurations. The dependence is very sim- 
ilar for the sites A, B and C, while S con f for the other 
two sites i.e. D and E depends on riFe in a way char- 
acteristic of the reference configuration. Eq. [3] shows 
that the configuration entropy, S CO nf, can be easily calcu- 
lated for any concentration of atoms on all sites, S CO nf = 

Sconf (XA ,X B ,X C ,X D ,Xe)- 



B. Free energy contribution 



The knowledge of AE, S mag and S con f for any 
concentration of Fe in the sublattices allows to spec- 
ify a formula describing a change of the free energy, 
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S m ag(xA,x B ,xc,x D ,XE)- Although changes AF = AE — T(S. 



in S mag due to changes in the individual site occupan- 
cies are monotonic, however, they deviate slightly from 



'mag + S con f) for any concentration 
of the alloying elements, and any temperature, AF = 
AF{xa ,xb,xc,xd,xe',T). Determination of the lowest 
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FIG. 3: (Online color) The configuration entropy, S CO nf, cal- 
culated versus the number of Fe atoms per unit cell, iiFe, for 
different atomic configurations of the er-FeV. The concentra- 
tion of vanadium is indicated, too. Scon/ -values determined 
for the occupation changes at each lattice site are plotted us- 
ing different colors. 



value of this function for the given temperature and con- 
centration, allowed to determine sublattice occupancies, 
Pi at the temperature T and Fe-concentration xp e — 
NiPi (where Nj denotes the occupation number 
for the sublattice i). The P^-values calculated for some 
chosen concentrations are shown in Fig. 4 for the temper- 
ature ranging between and 1800 K. At the first glace, 
one notices a significant site preference for the occupan- 
cies in the investigated phase. In particular, one can 
clearly see that the Fe occupancy of the sites A and D is 
actually complete and independent of the concentration 
and temperature. The occupancy variation with x and 
T is significantly larger for the other sites. It was found 
that increasing of x Fe in the alloy results in an increase 
of Fe content mainly on the sites C and E. For the latter, 
one can observe the highest sensitivity to temperature 
(Pe decreases with temperature) in comparison to much 
slower variations revealed for B and C site occupancies 
(that increase with temperature). 

It should be noted that the entropy contribution that 
predominantly effects the sublattice occupancies is the 
configuration entropy term. Neglecting S ma g in the cal- 
culation of Pi, depending on x and T as shown in Fig. 
4, leads to very small changes. It becomes clear when 
accounting for a similarity of all individual curves pre- 
sented in Fig. 2. There, S mag curves vary mainly with 
the total concentration of Fe/V atoms, and are very sim- 
ilar for all sublattices. This means that this form of the 
entropy does not substantially affect the site preferences. 
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FIG. 4: (Online color) The probability of finding Fe-atoms 
on the five lattice sites, Pi, versus temperature, T. The solid 
lines indicate the calculated values for different compositions 
of the cr-Feioo-z Vx with x ranging from 33 to 60 (the curves 
from top to bottom correspond to alloys with increasing a;). 
The vertical bars indicate experimentally found values of Pi 
in the same concentration range. 



The solid lines (Fig. 4) indicate the P^-values calcu- 
lated for different compositions of the u-Feioo-zV x with 
x ranging from 33 to 60 (full range of u-FeV occurrence). 
Vertical bars - marked in the plots - indicate values of Pi 
as determined at RT, but for the samples transformed to 
the (T-phase at 1300 K (that is why the bars are marked at 
this temperature in the plot). As the figure shows, the 
site occupancies calculated by the minimization of the 
free energy fairly agree with the values obtained from 
the neutron diffraction experiments. The agreement be- 
tween theory and experiment is less satisfying for the site 
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E, since the calculated P^-curves span over a wider range 
than the the width of the corresponding bar. 

The discrepancy observed between the calculated and 
experimental Pi-values is likely due to several factors. 
One of them is the simplified model used for the entropy 
calculations viz. without taking into account the phonon 
and the electron entropy terms. However, a possible in- 
fluence of these terms on final results seems to be rather 
small as previously found from the formation energy anal- 
ysis of the cr-phase in the Fe-Cr alloy system 5 . Unfortu- 
nately, we are not aware of any literature data reporting 
vibrational properties of the a— FeV that would allow to 
estimate the S p h contribution. On the other hand, bear- 
ing in mind a similarity of atomic and lattice parameters 
for Fe-V and Fe-Cr cr-phases, it seems reasonable to as- 
sume that the S p h contribution to the total entropy is also 
small in the present case. The effect of the S e i term was 
not taken into account mostly due to the fact that our 
discussion was focused on temperatures ~ 1300 K that 
are far above the Debye temperature of the investigated 
samples (~ 400 - 600 KM), 

Another factor that could be responsible for the dis- 
crepancy between the calculated and the experimental Pi 
- values is the assumed form (polynomials) of the func- 
tions for the AP and S ma g (hence, in consequence that of 
AP). Better fits were obtained by increasing the degree 
of the used polynomials, but in that case the discrepancy 
was not improved. 



IV. SUMMARY 

The formation energy, AP, magnetic, S magn , and con- 
figuration, S con f, entropy contributions were calculated 
and analyzed in the full range of the FeV-c-phase ex- 
istence using the electronic band structure calculations 
by means of the KKR method. It was found that AP- 
values strongly depend on the Fe concentration, and their 
variation observed for different site occupancies is charac- 
teristic of a given lattice site. The changes also strongly 
depend on the number of Fe atoms on the sites. Simi- 
lar, increasing dependence on the sample's composition 
exhibits the S magn term, but for a given composition it 
weakly depends on the lattice site. On the other hand, 
the composition dependence of the S con f term was found 
to be weak, although for a given composition, S con f, 
shows well-defined dependence on the lattice site, and for 
each site, a rather strong dependence on the number of 
Fe atoms occupying the site. The sublattice occupancies 
were determined for various compositions and temper- 
atures by minimizing the free energy. The occupancies 
of A and D sites were found to be almost independent 
of T and x, while the occupancies of the sites B, C and 
D showed a higher sensitivity to T and x. The KKR 
calculations combined with the analysis of the entropy 
contributions clearly showed that Fe atoms preferably 
occupy A and D sites in the FeV-er phase. This agrees 
well with the neutron diffraction measurements, and it 



confirms the general observation for a in binary alloy 
systems that the smaller atoms (in this case Fe) tend to 
mostly reside on the sites having the shortest distance to 
the nearest-neighbors (A and D), while the bigger atoms 
(here V) prefer to occupy the sites having more fee space 
to accommodate them (B, C, E). 
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Appendix A: Nomenclature 

There are several possibilities for description of sublat- 
tice occupancy in the FeV cr-phase. 

One can directly indicate rii (i = A through P), the 
number of Fe atoms occupying each site. Using that, 

n i — n Fe gives the total number of Fe atoms in the 
unit cell and n i gives the total concentration of iron. 

Sublattice occupancy by Fe atoms can be easily calcu- 
lated here 

i/riFe (hi other words Xi stands for 
the probability that randomly chosen Fe-atom belongs to 
the sublattice i, ^Xi = 1). 

On the other hand, P, = (Ni being the multiplicity 
of the i—th sublattice) denotes Fe-concentration on i-th 
sublattice (the probability, that randomly chosen atom 
on i-th sublattice is iron). 

Appendix B: Polynomial coefficients 

From the fitting procedure AP (in [eV/atom]), and 
Smag (in [eV/K/atom x 10~ 8 ] ) were determined, as fol- 
lows: 

AP = 0.17446 - 0.03178ni + 0.01269n 2 
-0.03404n 3 + 0.00768n 4 - 0.00689n 5 

+0.00258n^ 
+0.00226n 3 (ni + n 2 + n 4 + ra 5 ) 
-0.00160n 5 (m +7i 2 +«4 + n 5 ) (Bl) 



S mag = -1525.45 - 199.617m + 247.216n 2 
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+ 195.366n 3 + 624.922n 4 + 137.535n 5 
+62.3663n§ - 46.5917n^ - 11.4958^ 
-56.8196^2 + 117.04571? 
+287.426nin 2 + 35.9117nin 3 
-72.3616nin 4 + 0.05473nin 5 



+41.4487n 2 n 3 - 3.62722n 2 n 4 
-63.1291n 2 n 5 + 95.9298n 3 n 4 
-7.77130n 3 n 5 - 38.1694n 4 n 5 (B2) 
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